function y = func_w_for_integration(omega,chi,CONSTS)
% useless
    n_vec = (1:300)';
    m_vec = (0:20)';
    
    Wmn = zeros(size(m_vec,1),size(n_vec,1));
    Wm = zeros(size(m_vec,1),1);
    for wi = 1:size(omega,2)
        omega_i = omega(wi);
        for mi = 1:size(m_vec, 1)
            m = m_vec(mi);
            for ni = 1:size(n_vec, 1)
                n = n_vec(ni);

                if(chi)
                    Wmn(mi,ni) = func_Wmn_without_chi(omega_i,m,n,CONSTS);
                else
                    Wmn(mi,ni) =  func_Wmn(omega_i,1,m,n,CONSTS);
                end
            end        
            Wm(mi) = sum(Wmn(mi,:));    

        end
        y(wi) = 2*sum(Wm);
    end
    

end